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Abstract 

Results of lattice analysis indicate that the static potential in SU (3) gauge theory is 
proportional to eigenvalue of quadratic Casimir operator for the corresponding represen- 
tation with a good accuracy. We discuss the pattern of deviation from this Casimir scaling 
in gluodynamics in terms of correlators of path-ordered gauge-invariant operators defined 
on the worldsheet of the confining string. 



1 Introduction 



Physics of confinement in QCD has been attracting considerable attention since the Yang- 
Mills theory with dynamical quarks was established as a theory of strong interactions. Despite 
significant theoretical progress the operational framework for the exact calculations beyond 
perturbation theory has not been found, leaving aside numerical simulations of the theory on 
the lattice. In particular, one has no clear understanding of the relation between confinement 
and mass gap property. It is commonly believed that confining theory such as gluodynamics 
(i.e. QCD with no quarks) exhibits the mass gap, e.g. the lightest excitation over the physical 
vacuum is a massive particle. The latter one is indicated by the lattice calculations to be scalar 
++ glueball with the mass m ~ 1.5 Gev in the case of N = 3. On the other hand, if one studies 
confining potential between static sources in pure Yang-Mills theory, the physics of area law is 
governed by the formation of the confining string with nonzero tension \fa ~ 0.4 Gev and this 
process is not directly related to the details of glueball physics. The dynamical reasons for that 
are well known - corresponding effective interaction is suppressed by relatively large mass of the 
glueballs in the units of ^ l /a and also by 1/iV-factor. Put in simple terms, even if one acquired 
full knowledge of the glueball spectrum in gluodynamics, it would not help to compute generic 
Wilson loop. It is worth mentioning that in the abelian Higgs model with electrically charged 
condensate the mass of the photon (which plays the role of the lightest glueball in this theory) 
enters monopole-antimonopole string tension in rather nontrivial way and in combination with 
other quantities, giving the simplest example of the same nature. 

The standard approach of quantum field theory describes interactions in terms of par- 
ticle exchanges. In confining theories the applicability of this program is limited, in some 
sense: string creation, which is one of the most important manifestations of confinement can 
be understood neither in terms of gluon exchanges nor as a result of glueball interactions.^] 
Nevertheless it is desirable to establish some connection between different model approaches or 
at least to understand, where such relations are to be expected. The aim of the present paper 
is to make a step in this direction. It will be shown that the certain linear combination of the 
static potentials arising from the Wilson loops in different representations gets contributions 
from the exchange by colorless states. Our actual analysis will be presented for adjoint and 
fundamental loops which are of the main phenomenological interest but we believe that the 
mechanism is quite general. The paper is organized as follows: theoretical and lattice status 
of static potentials in different representations is briefly discussed in the Section 2, the Section 
3 presents the details and results of our analysis on more technical level, and short conclusion 
can be found in the Section 4. 

1 We do not discuss the models with the so called "confining gluon propagator" since they typically miss the 
phenomenon of gauge-invariant confining string formation as well. 
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2 Static potentials 



Recent lattice results on static potentials in different representations for SU(N) Yang-Mills 
theory fl], |], || [§] have attracted considerable interest (see 0] for references to early studies). 
It was argued in |5|, ||, |7[] that the data can serve as a good test for different phenomenological 
descriptions of confinement in gluodynamics and QCD. The reason for that lies in the impressive 
agreement between the lattice data and the so called Casimir scaling (CS) hypothesis M which 
states^] proportionality of representation-dependent confining potential (and, in more restricted 

(2) 

form, of the string tension) to an eigenvalue Co of quadratic Casimir operator C D for the given 
representation D. It is interesting that many microscopic models of confinement meet serious 
difficulties in reproducing CS (see discussion in 0). 

However, one is not to expect CS to be exact law in the case under consideration. As a 
matter of principle, CS should be violated at large distances because of screening. Qualitative 
picture for asymptotically large Wilson loops corresponds to the string tension gd equal to 
fundamental string tension a = Oj un ^ if the Wilson loop in the representation D nontrivially 
transforms under the center of the gauge group, or zero otherwise, in obvious contradiction with 
CS. On the other hand, the actual value of critical distance where effects of screening become 
important may be as large as 1.4 Fm J/J and one still has large enough "CS window", where 
the corrections to CS due to screening are exponentially small. As to perturbation theory, CS 
law for static potential is known to hold up to two loops in continuum theory [l(| however 
there is no any general theorem that CS should be exact property of perturbative potential 
and analysis of diagrams provides qualitative reasons in favor of CS violation at the three-loop 
level. It is also worth mentioning that CS is violated in the strong coupling expansion at the 
order 1/g 8 |l2j while it holds at lower orders. 



Another line of research confronts CS with alternative theoretical predictions for D-depen- 
dence of the static potential. The only law different from CS which has not yet been definitely 
ruled out in lattice simulations is the so called sin-formula, motivated by supersymmetric Yang- 
Mills theory analysis [13|, |j~4 |. Moreover, the authors of [15| claim that their numerical data 



for asymptotic string tension in d = 3 + 1 nonsupersymmetric SU(N) theory with N = 4, 6 



favor sin-formula and disagree with CS (see also |16|)- O n the contrary, authors of [17] claim 
their results to be in agreement with both CS and sin-formula predictions within two standard 
deviations in d — 3 + 1 Yang-Mills theory and definitely deviate from both in d — 2 + 1 
case. It should be noted that in all studied cases the numerical difference between CS law and 
sin-formula is actually quite small that plagues the discrimination between the two predictions. 

Unfortunately, as such, the sin-formula gives no nontrivial testable predictions for physi- 
cally interesting case of the gauge group SU(3). On the other hand, there are theoretical reasons 
to believe |T^] that it receives non-universal corrections. This seems to be a good motivation to 



take the predictions like CS or sin-formula as the first approximations and study the pattern 
of possible deviations from these regimes. Needless to say that despite numerically the two 



2 The term Casimir scaling was proposed in 
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predictions are surprisingly close, the qualitative physical picture behind them is drastically 
different. In any case, data from |L|, 0, |3|, |J as well as from |15|, [R| [17| undoubtedly indicate 



that CS as the first approximation for static potential works with percent-level accuracy at not 
too large distances. 

The starting point is the general expression for the Wilson loop 

(W D (C)) = (rr D P exp J A^dz^j ^ (1) 

where the generators Tg correspond to irreducible representation of dimension D. The nor- 
malized trace Tr D is defined as Tr D 1 D = j^Ti Id = 1, fundamental and adjoint generators are 
normalized according to 

TrT? und T b fund = 5 ab /2 ; Tr T^T b adj = N5 ab (2) 
The eigenvalue Co of quadratic Casimir operator Cp is defined as follows 

p{1) X rparpb rparpa ri 1 /q\ 

— °ab 1 D 1 D — 1 D 1 D — ^ D ' J-D 

We will confine our attention to adjoint representation in the present paper. The corresponding 
Casimir ratio d^/Dx = Cd 2 /Cd 1 is given by d a /f = 2N 2 /(N 2 — 1), where the fundamental 
Casimir Cf un d = (N 2 — 1)/2N and adjoint one C a dj = N equal to 4/3 and 3, respectively, for 
the case of SU(3). The static potential is formally defined as 

V D {/A,R) = -fin i hg(W D (C)) (4) 

where // stays for renormalization scale.^ The physical potential is given by 

Vg iys (R) = V D ( f i,R)-V D :elf ( f i) (5) 

where both terms in the r.h.s. diverge if \x —>■ oo while Vl, hys (R) is well defined in this limit (see, 
e.g. |33). It is known that perturbative series for the static potential in Yang-Mills theory does 



not exponentiate, contrary to abelian case [20]. Nevertheless it is convenient to exponentiate 



Wilson loop (P as (see review and references therein) 

(W D (C)) = exp ( £ i m AP(R,T)) (6) 



\m=l 



where each A^ 1 (R, T) is equal to well defined linear combinations of gauge-invariant field cor- 
relators (see (0), (0) below). We take (|]) as a nonperturbative definition of the potential in 
the present paper and problems of actual renormalization will not be considered, correspond- 
ingly the argument /i is omitted below. It is worth remembering however that all terms in the 
series defined by @ contribute to the potential and are to be renormalized accordingly. 

3 On the lattice the role of n is played by inverse lattice spacing a -1 
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In perturbation theory each A^f 1 (R, T) carries at least a factor g m which is small, the 
full series however is divergent due to factorial growth of the coefficients. On the other hand, 
nonperturbatively there are no general reasons to conclude that there exists some hierarchy be- 
tween correlators with different m. The general terminology calls the ensemble of the correlators 
stochastic if {A^ (R,T)\ > \A ( ™ +1) (R, T)\ and coherent, if \A { ™ ] (R,T)\ ~ |A^ +1) (i?, T)|. It 
can be argued that the vacuum of Yang-Mills theory is stochastic and not coherent [^, |23|, Ell 
and Casimir scaling is one of the strongest arguments in favor of that || |6|, [?|. The ultimate 
form of stochasticity is given by Gaussian dominance scenario when one keeps only the lowest 
nontrivial irreducible correlator in (|6]) corresponding to m — 2. The area law for the Wilson 
loop in the confining regime is given in the Gaussian vacuum in the following form: 

A { S{R,T)=d D/r af -R-T (7) 

where ~ (Tr F 2 (0))Tg and typical gluon correlation length T g is of the order of 0.2 Fm for 

pure N = 3 gluodynamics (see, e.g. [E_5], As it is clear from (0), Gaussian string tension 
is always proportional to the eigenvalue of quadratic Casimir operator for the representation D 
and therefore Gaussian dominance implies CS (not vice versa, however). Therefore to address 
the question of possible deviations from CS law, one has to go beyond bilocal approximation. 
This will be done in the next, somewhat more technical section. 



3 Deviations from CS law 



We are to study the general structure of CS violating terms in the cluster expansion. All 
calculations are performed for the simplest case of flat rectangular contour C = R x T in 
Euclidean space and its minimal surface. It does not mean any loss of generality since all 
expressions below can be trivially generalized to the case of arbitrary contour, whose minimal 
surface has the disk topology. In the case under study it is natural to adopt the so called 
contour gauge condition, introduced in [27, Ej| as generalization of Fock-Schwinger fixed point 
gauge f29| . It is defined by the condition that the phase factor orthogonal to the temporal axis 
is always equal to unity, i.e. 



$(x, y) — P exp 



i(y — x) / A(x(l — s) + ys, z±)ds 



= 1 



(8) 



where z 4 = x 4 = t/ 4 and A(x) = A a (x)T^. Without loss of generality one may also choose 



$(a; 4 ,y 4 ) = Pexp 



in 



i / A 4 (0, z 4 )dz 4 



X4 



= 1 



(9) 



The Wilson loop (|1]) takes especially simple form in this gauge 



(W D (C)} = Air D Pexp (i J A a 4 (z,t)T«dt 



(10) 
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where \z\ — R and the potential A\(z,t) is expressible in terms of the electric components of 
the field strength 



i 



Al(z,t) = J ' dszE a {sz;t) ■ E?(x) = F&(x) (11) 
o 

The contour gauges are well suited for studies of general properties of the gauge theories, but 
present considerable difficulties in perturbative calculations. Expanding (|H]) one gets, in the 
matrix form 



T tn-i \ 



(W D {C))=Tr D \i D + f^i n fdh... f dt n {M2,tn)-M*,h)) 

\ n=l I I 



(12) 



Each average in the r.h.s. of (112]) is proportional to the unit matrix in color space due to 
triviality of the vacuum quantum numbers, so the whole expression is, namely 

T tn-l 

Jdtx... J dt n (A A {z,t n )...A A {z,t x )) = U ■ (Tr D K%\R,T)) (13) 



With (H), (|) and ( p]) at hands one can easily rewrite every (Tr d (R> T)) in manifestly 
gauge-invariant form (see [21] and references therein). 

We will omit arguments R and T in the correlators (Tr o fC^(R, T)) below for simplicity 
of notation. We also do not discuss the actual dependence of each correlator on its arguments 
in the present paper. Lattice data suggest to fit the correlator as a sum of hard perturbative 
part and soft nonperturbative one, the latter being responsible for confinement. Notice that 
the series ( |T2"D is not strictly speaking perturbative series in the sense that each term in (p~2j) 
is given in perturbation theory by its own perturbative subseries. It is of no importance for 
our analysis which will proceed in terms of general expression and can be applied to whatever 
dynamics if it provides the set of correlators which obey Eq. ([13]). The perturbation theory 
itself can be regarded as a good example. 

It is now crucial to determine the relationship between (Tr o K,^) and A^K In the case 
of interest (Tr d^d) = 0' so the first term also vanishes while the relation between higher 
correlators is to be read off from the Taylor expansion of (0): 



Ag' 3) = (Tr D /Cg< 3) > ; Ag> = (Tr D /Cg } ) - - [(Tr D K%)\' (14) 
and the terms of higher orders are built in full analogy with ( |T3D : 

A^ = (Tr D /4 m) > -(...) (15) 
where (...) stays for the sum of products of the correlators (Tr d^d ) °f lower orders, p < m. 



1 



-(2)\ 



2 



The expression (15) does not represent the Green's function of any physical colorless state, 
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instead, it describes the propagation of some gluelump state, which is not a physical excitation 
in gluodynamics. It is also worth saying that (|I4]) , (|15|) should be considered as a definition, 
valid beyond perturbation theory (each correlator is given by the corresponding perturbative 
series in the latter). 

The potential Vd(R) is given up to an additive constant by 

1 oo 

Vb(J2) = - T Km - £z m A^ (16) 

m=2 

We are going to study now the structure of special linear combination of potentials, measuring 
the deviations from Casimir scaling 

SV D2 / Dl = V D2 (R) - d D2/Dl ■ V Dl (R) (17) 

We confine our attention to the case of fundamental and adjoint representations in this paper, 
in this particular case ( |T7j ) takes the form 

-i oo / 9 at2 \ i oo 

SV a/t = - Hm i £ f A« - ^ ■ A" = - Um I £ roA$ (18) 
As a next step, we are exploiting the identity 

<l^)| 2 ) = ^+(l-^)(^)) (19) 

We need to compute the quantity (jWf(C)\ 2 "j in terms of correlators flTB|). It can be done as 
follows. The average of the product is given by 




where the trace Tr 12 = Tr / x Tr f 2 acts on the color indices corresponding to the first and 
the second loop independently. Global gauge invariance dictates the following form of the 
corresponding averages (see also p0|): 

()Cf) = i r (Tr f ICf) ; (K$ K%) = 
+ e ■ j£- ((Tr , K$) - (Tr Tr ,X#>) (21) 
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The matrices 1/, 1 and e can be written in index notation as 

[l/]ai/3i = ^ail/31 [l]ai/3i;a 2 /32 = <W /3i ^0:2 ft> > [^\a 1 /3 1 ;a 2 P2 = ^ai/3 2 ^a2/3i (22) 

where ai i2 ,/3i j2 are fundamental indices running from 1 to N. The algebra they obey encodes 
the effects of path ordering. In our case since the orientations of minimal surfaces corresponding 
to the first and the second multiplier in ( ^0|) are antiparallel, the matrix e should be multiplied 
from the right with respect to indices carrying subscript "2" but from the left with respect to 
indices carrying subscript " 1" . It results in the following algebra: 

g .1 = i . g = e ; P = i; e 2 = Ne (23) 

and the traces are given by 

Tr 12 i = l ; Tr 12 e=l (24) 
It is convenient to represent (p0| ) in the following form: 

(w f (C)Wj(C)) = Tr 12 exp (i ■ A + e • B) = Tr 12 (i ■ exp A+ 

+ ± . ( exp ( A + BN) - exp A) J = ^ exp(A + BN) + (l - -^j exp A (25) 

Comparing (|25]) with fll9|), (|20| ) and ( pTj) one concludes that 

exp(A + BN) = 1 

oo oo 

exp A = (W a (C)) = 1 + E * n < Tr /4 n) ) + E("0 P (Tr //CjV 

n=l p=l 

oo oo -n/ -\p 

+ EE U^T (Tr Tr > - < Tr / 4? 4?)) (26) 

n=lp=l iV 1 

This is exact expression for adjoint Wilson loop average, written in terms of fundamental 
correlators and traces and as such it is suitable for our analysis. We are to study the difference 

5V a/f = -fim i (log(W a (C)) - J^-log(W f (C))] (27) 



where the first term in the curly brackets in fl2"7|) is equal to log(W a (C)) = A where A is defined 
by (p6|), while the second term is given by formal series 

k 



Coo \ oo / I \k 

l + E^(Tr/4 n) ) =-ES^ 
n=l / fc=l K 



£^<Tr/4 n) > 
n=l 



(2f 



To proceed further one can exploit the identity ([W f (C)] a p[W}(C)]p a ) = 1 or, in terms of 



correlators 



E *"<Tr fICf) + £H) P (Tr ) = -£E T(-^(Tr / 4?) (29) 

n=l p=l n=l p=l 
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It is now straightforward to rewrite ( ^6|) as: 

N 2 / oo 
i«=l 



log(W a (C)) = A = log I 1 + I £ *"<Tr + E(-0 P (Tr //Cf ) + 



-(p)\ 



p=i 



-E 



oo m— 1 

+ EE i i- i )^/4r ) ^/^ ) > 

m=2 p=l 



;-i) fc / at 



2 \ fc fc 



\ iV 2 - , 
fc=i ^ V v x / 2=0 



E^m/4 nj ) + E(-^/^f> 
p=i 



n=l 



oo m— 1 

53 r ^ (-l)P(Tr ; /C^- p) Tr } Kf) 

m=2 p=l 



k-l 



(30) 



Till now we have not made any approximations yet and the expressions (^), ( |3"0"D are valid 
for arbitrary contour C and any contour gauge. It is worth reminding that despite each term 
in (1280 , (|30|) is gauge- invariant, its actual value depends on the profile of the contours used in 
(i)> @ but this dependence is cancelled after summation of all terms. Likewise the property of 
path-dependence of Gaussian correlator (found in to be rather strong) shed no light on the 
problem of CS since gauge- invariant path-dependent bilocal correlator (Tr D F^ u (x, x )F pa (y, x )) 
where F^x, x ) = ^(xQ,x)F fMU (x)[^(x ,x))^ and $(x ,x) are the phase factors (chosen to be 
unit matrices in our analysis) is always proportional to Co for any path, taken in $(xo,x). 

We can now use the Euclidean rotational 4-invariance of our problem which dictates 
(W f (C)) = (W}(C)) or, at the level of correlators, 



^^(Tr / /C^) = 5:H) P (Tr / /Cf) 
P =i 



(31) 



71=1 



Let us now look at the general structure of terms in 5V a /f. It is straightforward to notice 
V/ 

-(m-p) rr, y{p)\ 



that each 5A^l gets contributions from irreducible correlators of gauge-invariant operators of 



the type (Tr fK,\ Tr jK-JJ). The contribution from the single correlator of the same order 

can be shown to be exactly cancelled in ([27D (this cancellation happens in the terms 
corresponding to k = 1 in fl28|) and (0)). To illustrate the point, consider the lowest nontrivial 
order of cluster expansion violating CS, which corresponds to m = 4. One gets 



N 2 



N 2 - 1 



1 + 



N 2 - 1 



(Tr//Cf) 



(32) 



This expression is gauge-invariant and the fields A±(x) in the correlators (Tr/C^) are functionals 
of the field strengths F i4: (x) as given by (|TT|). 

The l.h.s. of expression (^) vanishes identically if CS law is exact. Posed differently, it 
measures the deviation from CS in terms of integral moments of particular correlators, entering 
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the r.h.s. of fl32]) . The coefficient in front of the second term in the r.h.s. of (^) (equal to 5/4 
in the case of SU(3)) provides exact cancellation of reducible contributions from the first term, 
so both sides of equality (^) have area law asymptotics in confinement regime. In large N limit 
they both vanish, as it should be. The same pattern holds in higher orders of cluster expansion. 
Casimir scaling violating contributions are given at fixed m by the irreducible correlators of 
colorless operators (Tr //Cy™ ^ Tr fK$)- ^ * s ^° ^ e compared with the expressions flT4|), (P~5|) 
where the potential itself is written in terms of correlators (Tr fKY 1 ). 



4 Conclusion 



Expression (^2|) and analogous formulas for higher orders of cluster expansion formalize the 
statement made in that violations of CS could be a result of exchange by colorless states. 
We have considered only the case of adjoint and fundamental potentials in the present paper 
but no reasons for breaking of this mechanism in other representations are seen. The qualitative 
picture corresponds to the confining string worldsheet populated by irreducible correlators like 
(Tr £) k!§) and correlators of higher orders made of gauge- invariant expectation values of path- 
dependent field strength operators. On the minimal surface the whole ensemble contributes to 
the total string tension in such a way that Gaussian contribution is numerically dominant, as 
CS suggests. In other words, CS indicates that this ensemble is quasi-free in the sense that 
interactions between such gauge-invariant objects as Tr FF mediated by glueball states and 
destroying CS are rather weak. Phenomenologically, we expect the r.h.s. of expressions like ( p2l) 
to be suppressed by some power of inverse glueball mass. The actual smallness of deviations 
from CS as observed on the lattice simply follows in this case from the fact that even the lightest 
glueball is rather heavy, of the order of 1.5 GeV. Moreover, one can have a glimpse of general 
reason why many microscopic confinement models fail to reproduce CS. To this end, the model, 
apart from explaining string tension formation and area law of the Wilson loop at intermediate 
distances for all representations (" CS plateau" ) must also describe the phenomenon of mass gap 
itself and come out with the correct (rather large) masses of colorless states. One can mention, 
on the other hand, that reasonable prediction of the instanton liquid model for the mass of the 
lowest glueball (1.4 ± 0.2) GeV |52|| is in some contradiction with the strong violations of CS 



for the static potential in this model . This paradox is possibly resolved by the fact that the 
instanton-antiinstanton ensemble is coherent and the cluster expansion is not well suited for 
analysis in this case. 
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